Nature of carriers in one-dimensional correlated electron systems 
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Electron densities at the edge site in one-dimensional Hubbard model and t-J model are studied 
by using the Bethe ansatz solutions and exact diagonalization method. It is found that the boundary 
is electron-attractive, or equivalently hole-repulsive, near half-filling. We propose a new criterion to 
determine whether the carriers in a strongly correlated system are electron-like or hole-like depending 
on the electronic behavior at the boundary. 



A suitably substituted impurity atom in condensed 
matter systems can play a role as an atomic scale probe 
for local response of the system. In particular in strongly 
correlated electron systems, the response to an impurity 
offers a key to understand the nontrivial properties of 
the system. Therefore, a large number of studies have 
been carried out on the impurity effects in the strongly 
correlated electron systems, especially in high- Tc super- 
conductors. One of the most important findings concern- 
ing the impurity effects in high- Tc superconductors is the 
appearance of local magnetic moments around a non- 
magnetic impurity, such as Zn or Li atom, substituted 
for Cu in the Cu02 plane in underdoped compounds, 
diiii. 

Recently, we have pointed out that whether the impu- 
rity in the Cu02 plane tends to expel electrons or holes 
is crucial for the appearance of induced moments or the 
local enhancement of antiferromagnetic (AF) correlation 
. If the system is in the overdoped region, it has been 
shown that electrons are expelled from the neighborhood 
of the impurity and thus a hole-rich region is locally 
formed around the impurity. As a result, AF correla- 
tion is collapsed because the effective exchange coupling 
between electrons is locally suppressed. Conversely, if 
the nature of the carriers is hole-like in the underdoped 
region or close to half-filling, we expect that the holes are 
scattered by the impurity and then an electron-rich re- 
gion is formed around the impurity. If this is the case, it is 
natural to consider that local magnetic moments appear 
near the impurity since the local region approaches half- 
filling. Therefore it is very important to clarify the nature 
of carriers in correlated electron systems, on which the 
local magnetism around an impurity strongly depends. 

In the two-dimensional systems, however, it is difficult 
to determine the nature of carriers doped in the Mott in- 
sulator. Slave-boson mean-field theory suggests that the 
charge degrees of freedom can be represented as a small 
number of holons, however this naive picture is uncer- 
tain because there exists the large Gauge fluctuation. In 
this paper, we investigate the nature of carriers in one- 
dimensional (ID) systems where the essence of strong 
correlation can be clarified in exact solutions. 



It is worthwhile noting here that even in ID systems, 
the nature of carriers has not been understood because 
there is no Hall effect. Thermopower was suggested to 
clarify the charge of the carriers, but its calculation is 
only phenomenological 0- Instead of these quantities 
we propose that the local electron density near the im- 
purity gives a good quantity for determining whether the 
carriers are hole-like or electron-like. When the impurity 
is repulsive, it scatters the mobile carrier. Near half- 
filling, we expect that the holes are scattered and the 
hole density should decrease around the impurity. Using 
this criterion we determine the phase diagram in the ID 
Hubbard model. 

When we assume the impurity scattering is in the uni- 
tarity limit, the impurity is equivalent to the bound- 
ary in ID systems. We calculate the electron den- 
sity at the edge site of the ID Hubbard model and 
t-J model with boundar y by using the Bethe ansatz 
solutions 0, 0, 0, ^2 Q| and the exact diagonal- 
ization method. 

The Hamiltonian of the open Hubbard chain with 
boundary fields, pi and pl, is given as 
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following the notation in ref. 9]. The hopping integral, 
has been taken as a unit of energy {t = 1). We describe 
the length of the chain by L, which is assumed to be 
an even integer. The average number of electrons at the 
edge site (ni) taken with respect to the ground state |<I>o) 
is obtained from 
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The Bethe ansatz equations for the above Hamiltonian 
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where u = U/A, N is the total number of electrons and 
M is the number of down spins In this letter, we 
only discuss the repulsive Hubbard model. For this case 
it is useful to denote the negative fcj 's and A^'s as fc_j = 
—kj, A_/3 = —A/3 . Using standard procedures, the Bethe 
ansatz equations for the ground state can be rewritten 
as the following linear integral equations for (T£(fc) and 
tT£(A), which represent the densities of solutions for the 
real quasimomenta kj and spin rapidities A/3: 
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where 



are the standard kernel for the Hubbard model, and 
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The values of fc+ and A"*" are determined from the condi- 
tions 
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In one-dimensional systems with open boundaries, the 
first finite-size correction for the ground state energy El 
is of order L*^, which is the contributions from the bound- 
ary. Thus the ground state energy of the present model 
can be written as 



El = Le^{k+,X+) + ei(fc+, A+) + 0(1/L) (7) 



with 
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Here, edk) and £s(A) are the dressed energies defined as 
£c{k) = —2cosk + fj, 

+ — / dA'i<ri(sin/c- A')e.(A), (9) 
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If the energy El in the thermodynamic limit is minimized 
for fc+ ~ ko and A+ = Aq, the parameters ko and Ao are 
obtained by 
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These conditions are equivalent to 

Seiko) = 0, £s(Ao) = 0, 
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which are the same as in the periodic-boundary case. 

Now we expand Coo and ei up to the second order in 
Afc = fc+ — fco and A A = A"*" — Ao. As far as the conditions 
eq. Hll|) are satisfied, we find Afc and AA are of the order 
of 1/L, and fc+ and A''" in eq. Q can be replaced by fco 
and Aq, whose values are determined from the conditions 
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Here cr^(A;) and cr^(A) are obtained by taking the limit 
of L ^ oo in eqs. I@J and Since Aq = oo for the case 
of = 2Af, (7^ (A) can be solved from the Fourier trans- 
formation. Substituting cr^(A) into the integral equation 
for cr^(fc), we obtain 
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Finally to obtain ei(fco,Ao), we need £c{k) and es(A). 
In the same way as (T§^(fc), we have 
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As a result, the electron density at the edge site is given 
by 
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In the actual calculation, we solve eqs. (|13|) and H14|) to 
obtain ko and crj^(fc) for each value of U/t and N/L. 
Then using the obtained feg, we solve Sc{k) in eq. lfTB|l to 
calculate ni. 

Figure 1 shows the obtained ni as a function of the 
electron density, n — N/L, of the whole system. In the 
noninteracting case {U/t = 0), ni is always smaller than 
n in less-than-half- filling case {n < 1), which means that 
the carriers are electrons. For n > 1, tt-i can be easily 
obtained from the electron-hole symmetry and is larger 
than n, indicating that the carriers are hole-like. As U/t 
increases, the crossover point from ni > ntoni < n 
approaches n = 0.5 from n — I. If we take the limit of 
U/t —>■ oo, the ni-n relation is just half-size of that for 
U/t — 0, i.e., n = 2 for U/t — corresponds to n = 1 
for U/t — > oo. This can be easily understood because 
the charge degrees of freedom are expressed as a Slater 
determinant of spinless fermions |l4| . Thus, n — 0.5 
means the half- filled band of the spinless fermions. 

Figure 2 shows the crossover points as a function of 
t/U. Apparently, in the region of n 1 and large U/t, 
nil is larger than n. From these results, we interpret 
that the holes are scattered by the impurity (in this case, 
by the boundary) and thus the nature of carriers in this 
region is hole-like. This is one of the typical effects of 
strong correlation in the doped Mott insulator. 

Next we study the same problem in the t-J model. The 
Hamiltonian of the open t-J chain with boundary fields 
is given as 
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FIG. 1: The average number of electrons at the edge site of 
the open Hubbard chain as a function of the electron density 
of the whole system. A dashed line indicates n\ = n line. 
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FIG. 2: The phase diagram in ID Hubbard model, 
nature of carriers is hole-like in the shaded area. 
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following the notation in ref. [l^- The Bethe ansatz 
equations for the open t-J model are available only when 
J /t = 2 under periodic 0, and open boundary con- 
ditions 0, 0| . Following the same procedure as for the 
Hubbard model, we obtain the electron density at the 
edge site. For the other values of J/t, we calculate ni in 
the exact diagonalization of the system with 16 sites. 

Figure 3 shows the results of the t-J model in the same 
manner as in Fig. 1. Note that J/t = corresponds to 
U/t ^ oo in the Hubbard model. The exact diagonal- 
ization (ED) results show that, as J/t increases, ni de- 
creases for any values of n and hole-like region shrinks. 
This behavior is consistent with that of the Hubbard 
model in the large U region. For J/t = 2, the ED re- 
sults slightly overestimate the value of ni compared to 
the Bethe ansatz (BA) results. This is the finite size ef- 
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feet, i.e., there is no enough space where electrons escape 
from the boundary in small systems. 

It is interesting to compare the results for the t-J 
model with that for the Hubbard model. For J/t = 2, the 
BA result shows that ni < n for any values of n, which 
means that the carrier is always electron-like. Actually, 
the result for J/t = 2 is very close to that for the free 
fermion case {U/t = in Fig. 1). The reason is as fol- 
lows. When J/t = 2, the electron motions driven by the 
hopping term and by the exchange term have the same 
amplitude t = J/2. Therefore, electrons can move almost 
independently with each other and thus a free electron- 
like behavior can be seen for J/t = 2 [13. Furthermore 
the result for J/t = 1 is close to that for U/t = 4, which 
is consistent with the relation between J and [/, i.e., 
J = 4t^/[/. Even though the so-called three-site terms, 
which appear in the perturbation theory with respect to 
t/U , are neglected in eq. H18|l . we can see that the na- 
ture of carriers has a strong correspondence between the 
Hubbard model and the t-J model when J/f^ 1. 
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FIG. 3: The average number of electrons at the edge site 
of the t-J model as a function of the electron density of the 
whole system. The results are obtained by using the exact 
diagonalization (ED) or the Bethe ansatz equations (BA). A 
dashed line indicates n\ = n line. 

In summary, we have investigated the electron number 
at the edge site of the one-dimensional Hubbard model 
and t-J model by using the Bethe ansatz solutions and 
the exact diagonalization method. We propose a new cri- 
terion to determine whether the carriers in strong corre- 
lated systems are electron-like or hole-like, which is usu- 
ally difficult in ID because of the lack of the Hall effect. 
We expect the similar behavior in 2D Hubbard model or 
t-J model. Although reliable estimation of the crossover 



point will be difficult, it is natural to expect that the 
carriers are hole-like (or holons) near half-filling, which 
causes the local magnetic moments around impurities. 
As the system approaches half-filling, the number of the 
hole-like carriers decreases leading to the metal-insulator 
transition. 
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